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By removing one empty site between two occupied sites, we map the ground states of chains of hardcore 
bosons and spinless fermions with infinite nearest-neighbor repulsion to ground states of chains of hardcore 
bosons and spinless fermions without nearest-neighbor repulsion respectively, and ultimately in terms of the 
one-dimensional Fermi sea. We then introduce the intervening-particle expansion, where we write correlation 
functions in such ground states as a systematic sum over conditional expectations, each of which can be ul- 
timately mapped to a one-dimensional Fermi-sea expectation. Various ground-state correlation functions are 
calculated for the bosonic and fermionic chains with infinite nearest-neighbor repulsion, as well as for a ladder 
model of spinless fermions with infinite nearest-neighbor repulsion and correlated hopping in three limiting 
cases. We find that the decay of these correlation functions are governed by surprising power-law exponents. 

PACS numbers: 71.10.-w, 71.10.Fd 



I. INTRODUCTION 



Exact solutions hold a special place in the theoretical con- 
densed matter physics of interacting electron systems. Al- 
though they can be obtained only for very specific models, 
these proved to be very useful in understanding the behaviour 
of more general models of interacting electrons, or informing 
us of novel physics that we would otherwise not suspect from 
approximate treatments. In particular, our present paradigm 
of two universality classes, Fermi liquids versus Luttinger 
liquids, for low-dimensional systems of interacting fermions 
came out of exact solutions showing separation of the charge 
and spin degrees of freedom. 1,Z ' 3 ' 4 

In this paper, we report further surprises coming out of 
the exact solution of models of hardcore bosons and spinless 
fermions with infinite nearest-neighbor repulsion^ We con- 
sider chain models 
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of hardcore bosons (U — > oo) and spinless fermions, as well 



as a ladder model 
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of spinless fermions. In this ladder model,-- the correlated hop- 
ping —t'Cj jNi>&,j+iCij + 2 is the simplest term we can introduce 
to blatantly favor the emergence of superconducting order. 

Throughout this paper, we will specialize to the limit of in- 
finite onsite repulsion U — > oo and infinite nearest-neighbor 
repulsion V — > oo. More precisely, we admit only config- 
urations in which each site can be occupied by at most one 
particle, with no simultaneous occupation of nearest-neighbor 
sites. We will show how the nearest-neighbor excluded chain 
models can be mapped to the nearest-neighbor included chain 
models 

H%*> = -t£ [b)b j+i + b) +1 bj\ + uY j n j (nj - 1), 

j 

and ultimately solve for the ground states of the former in 
terms of the one-dimensional Fermi sea. We will also show 
how the ladder model can be solved exactly in three limiting 
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cases, by mapping their ground states to those of the chain 
models given in Eq. ( 11.11 ) and Eq. ( 11.31 ). These analytical re- 
sults were used to guide a density-matrix analysis of correla- 
tions for the ladder model, first using the exactly diagonalized 
ground states,? and later using the density-matrix renormal- 
ization group.— 

For the rest of this paper, we will consistently use uppercase 
letters Bj and B\ (Cj and Cj) to denote hardcore boson (spin- 
less fermion) annihilation and creation operators on nearest- 
neighbor excluded chains, and lowercase letters bj and b\ 
(cj and c\) to denote hardcore boson (spinless fermion) an- 
nihilation and creation operators. Similarly, Nj = B^Bj (or 

Nj = C\Cj) and rij = b^bj (or nj = etc/) are the hard- 
core boson (spinless fermion) occupation number operator on 
the nearest-neighbor excluded and nearest-neighbor included 
chains respectively. Hereafter, we will also use excluded to 
refer to all quantities associated with the nearest-neighbor ex- 
cluded chain, and ordinary to refer to all quantities associated 
with the nearest-neighbor included chain. 

Our paper will be organized as follows: in Sec. [Til we 
will describe an analytical map that establishes a one-to- 
one correspondence between the Hamiltonian matrices of 
the excluded and ordinary chains of hardcore bosons and 
spinless fermions, before developing a systematic expansion 
that would allow us to calculate ground-state expectations in 
bosonic and fermionic excluded chains. We then present and 
analyze in Sec. [Ill] correlation functions calculated using the 
analytical tools developed in Sec. [n] for excluded chains of 
hardcore bosons and spinless fermions. Following this, we 
write down in Sec. [IV] the exact ground states of the ladder 
model given in Eq. ( 11.2b in three limiting cases, and calculate 
various ground-state correlation functions, before summariz- 
ing our results and discuss the interesting physics they imply 
in Sec.[V] 

II. MAPPINGS AND TECHNIQUES 

In Sec. Ill Al we establish a one-to-one correspondence be- 
tween states of the nearest-neighbor excluded and nearest- 
neighbor included chains. We explain how the Hamiltonian 
matrices of the two chains, and hence their energy spectra, are 
identical to one another. In the infinite-chain limit, we then 
show how we can write the ground state of the excluded chain 
in terms of the ground state of the ordinary chain, and ulti- 
mately be written in terms of the one-dimensional Fermi sea. 
In Sec. Ill Bl we show how the ground-state expectation be- 
tween two local operators can be calculated for the excluded 
chain, by writing it as a systematic sum over conditional ex- 
pectations, each of which associated with a fixed configuration 
of intervening particles. 

For the sake of definiteness, let us consider open chains of 
a finite length L and total particle number P. Sites on these 
chains are indexed by j = 1 , . . . , L. Since the models in 
Eq. ( 11.11 ) and Eq. ( 11.31 ) conserve P, the infinite-chain limit is 
obtained by letting L — » oo, keeping the density of particles 
N - P/L fixed. Ultimately, the results we present in this sec- 



tion will not depend on what boundary conditions we impose 
on the chain (which is what we would expect in the infinite 
chain limit). 

For convenience, we establish some notations to cover bo- 
son and fermion cases together in the same formula. Let us 
call 

\J)=A)A\ 2 ---A) p \0) l , (2.1) 

an excluded configuration, where A — B for hardcore bosons, 
A = C for spinless fermions, and the sites < j\ < ji < 

■ ■ ■ < jp < L are such that j p+ \ > j p + 1. We will also employ 
the labels a and B for distinct P-particle configurations |J") 
and |J^), i.e. the P-particle configurations {/*, j", . . ., j a p ) and 

{f v f 2 ,...,f p ) differ in at least one site. 
Similarly, let us call 

*>-«]«] •••«;. <>>,. (2.2) 

an ordinary configuration, where a - b for hardcore bosons, 
a — c for spinless fermions, and the the sites < j\ < ji < 

■ ■ ■ < jp < L are such that j p+ \ > j p + 1. The labels a and 
B will again denote distinct P-particle configurations and 

We will also consistently denote the Hamiltonian of an 
excluded chain by Ha, where Ha - H£jv for hardcore bosons, 

(c f) 

Ha = H)y for spinless fermions, and the Hamiltonian of an 
ordinary chain by H a , where H a = fljff for hardcore bosons, 

(c f) 

H a - H) ' J for spinless fermions. 

A. Mapping Between the Excluded and Ordinary Chains 

In this subsection, our goal is to establish the one-to-one 
correspondence between states of the excluded and ordinary 
chains, and to show that as matrices, the Hamiltonians dl.lb 
and (II. 31 are identical. To do this, let us note that an excluded 
chain with L sites has fewer P-particle states than an ordinary 
chain of L sites, because of the infinite nearest-neighbor re- 
pulsion. Therefore, we can form a one-to-one correspondence 
between excluded and ordinary states only if the length U of 
the ordinary chain is shorter than L. There are several ways to 
systematically map excluded configurations to ordinary con- 
figurations: we can (i) delete the site to the right of every par- 
ticle, if it is not the rightmost particle; or (ii) delete the site to 
the left of every particle, if it is not the leftmost particle. We 
can easily check that these maps produce the same ordinary 
configurations for finite open chains. We expect this to hold 
true even as we go to the infinite chain limit. For the rest of 
this paper, we will adopt right-exclusion map 

'U-Vr-V/' i o>^, v ,;....„ ;; o> (2.3) 

that maps a P-particle configuration on an open excluded 
chain of length L to a P-particle configuration on an open or- 
dinary chain of length L' . The empty site to the right of each 
occupied site in the open excluded chain is deleted, to give a 
corresponding configuration for an open ordinary chain^ As 
illustrated in Fig. [T] we do not delete any empty site to the 
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right of the Pth particle, and thus the effective length of the 
open ordinary chain is L' = L - P + 1. This nearest-neighbor 
exclusion map was first used by Fendley to map a supersym- 
metric chain of spinless fermions to the XXZ chain. 9 It tells 
us that an excluded chain with density 



_ P 
N =L 



gets mapped to an ordinary chain with density 



N 



U L-P+l 1-N + (1/P) 
Thus, in the limit of L, P — > oo, 

N 



l-N 



(2.4) 



(2.5) 



(2.6) 



FIG. 1: Schematic diagram illustrating how we map from P-particle 
configurations on an open excluded chain of length L to P-particle 
configurations on an open ordinary chain of length V = L — P + 1, 
by deleting one empty site to the right of a particle that is not the 
rightmost particle. 

For P-particle excluded configurations \J a ) and |J^), the 
matrix element (J"\Ha\J^) is nonzero only when |J a ) and U' 3 ) 
can be obtained from one another by a single particle hop- 
ping to the left or the right. When this is so, the ordinary 
P-particle configurations \j a ) and they map to under the 
right-exclusion map are also related to each other by a single 
particle hop. Thus, we have 



(3 a \H A \f) = -t = Q a \H a f). 



(2.7) 



Since there is a one-to-one correspondence between P-particle 
configurations on an open excluded chain and P-particle con- 
figurations on an open ordinary chain, Eq. ( 12.71 ) tells us that 
Ha and H a are identical as matrices in their respective con- 
figurational bases. We therefore conclude that the P-particle 
energy spectra of the two chains coincide, and that there is 
a one-to-one correspondence between the energy eigenstates, 
ll^f)} for the excluded chain, and {]¥')} for the ordinary chain. 
That is, if |J) i — > j), and \*¥) h* then |J> and j> have 
the same amplitudes in |*f) and |T") respectively. This result 
has profound implications on the thermodynamics of the two 
chains, as well as that of the ladder model in the three limiting 
cases described in Sec. IIV Al because their partition functions 
are the same. However, for the rest of this paper, we limit 
ourselves to the ground states of the infinite excluded and or- 
dinary chains, as well as those of the infinite ladder. 



B. Ground-State Expectations of the Excluded Chain 

In this subsection, we explain how the expectation {O) of 
an observable O in the ground state of the excluded chain 
can be computed, by relating it to the expectation (O') in the 
ground state of the ordinary chain, for an appropriately chosen 
observable O' satisfying some basic correspondence require- 
ments that we shall outline. Specifically, we are interested 
in the correlations {0\Ot) between two local observables 0\ 
and Oo separated by a distance r within the excluded-chain 
ground state. However, the right-exclusion map maps ex- 
cluded matrix elements (JHOiC^U^) to ordinary matrix el- 
ements (S a \0\0' 2 \f) in which 0\ and 0' 2 are separated by 
varying separations. To deal with this problem, we develop 
a method of intervening-particle expansion involving a sum 
over conditional expectations. 

To begin, let us look at the ground states 



PP > = 2>o(J)IJ>, 



I.] I 



m> = 2>oCi)Li> 

U! 



(2.8) 



of the excluded and ordinary chains respectively. To take 
advantage of the equality of amplitudes, *Fo(J) = ^oG) if 
U) lj) under the right-exclusion map, we want (j^lO'lj 8 ) 
to have a simple relation with (J a \0\jP). While it is possible 
to pick a such that {j a \0'\f ) = <J"|<9|J^> for all a and /?, we 
find that it is more convenient to pick O' such that 



N n 



(2.9) 



For example, if O — Nj = CjCj, we can pick the correspond- 
ing observable to be O' = nj = CjCj, in which case we find 
that 



{3 a \Nj\f)=NS afi , 



(2.10) 



which satisfies Eq. ( 12.91 ) trivially. We call O and O' a corre- 
sponding pair of observables, if Eq. (12.91 ) is satisfied for all a 
and allowing us to write the very simple relation 



<o>= l -{0') 



(2.11) 



between their ground-state expectations. 

Since we are mostly interested in correlation functions 
within the excluded chain ground state, let us look at expec- 
tations of the product form (OjOj +r ), where Oj acts locally 
about site j, and Oj +r acts locally about site j + r. Becasue 
the number of particles p between sites j and j + r varies from 
excluded configuration to excluded configuration, these sites 
get mapped by the right-exclusion map to sites on the ordi- 
nary chain with varying separations r — p. Therefore, to cal- 
culate the excluded chain ground-state expectation (OjOj +r ) 
in terms of ordinary chain ground-state expectations, we first 
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write down an intervening-particle expansion 

(OjO j+r ) = J](OjQ p Oj +r }, (2.12) 
p 

where (OjQ p Oj +r ) are conditional expectations. Here p is 
a vector of occupation numbers within the intervening sites, 
and Q p is a string of factors, each of which is either Nj+ S or 
(1 - Nj+s), 1 < s < r — 1. The sum is over all possible ways to 
have intervening particles between Oj and 0j +r . For each ex- 
cluded term (OjQ p Oj +r ) in Eq. ( 12.121 ), we then write down the 
corresponding ordinary expectation {0'-Q' v ,0'j +r p }, and there- 
after sum over all corresponding ordinary expectations, 

(OjO j+r ) = jJ](0' j Q' 1 ,0' j+r _ p ), (2.13) 
p' 

making use of Eq. ( 12.1 11 1. The vector p' of occupation num- 
bers is obtained from p using the right-exclusion map, and 
contains the same number p of occupied intervening sites. 

To illustrate how the corresponding expectations 
(0'jQ' p ,0'j +r p ) can be constructed, let us write Eq. d2.121 > out 
explicitly as 

(OjO j+r ) = (Oj(t - N j+1 )- • • (1 - iV i+r _i)0 7+r > + 
(O j N j+l ---(t-Nj +r - l )Oj +r ) + --- + 
(Oj(l-N J+1 )---N j+r -iOj +r ) + 
(OjNj +l Nj +2 ■■■(!- N j+r -i)O j+r ) + ■■■ + 
(Oj(t - N j+l ) ■ ■ ■ N j+r - 2 Nj +r -iO j+r ) + ■■■ + 
(OjNj +l Nj + 2---Nj +r - l O j+r }, 

(2.14) 

each of which contains intervening particles at fixed sites. 
We call terms in the expansion with p intervening N/s 
the p-intervening-particle expectations. Because of nearest- 
neighbor exclusion, most of the terms in Eq. ( 12.141 ) vanish. 

Next, we map each conditional excluded expectation in 
Eq. (12.141 ) to a corresponding conditional ordinary expecta- 
tion following the simple rules given below: 

1 . Nearest-neighbor exclusion. To ensure that we do not 
violate nearest-neighbor exclusion, we make the assign- 
ment 



A , A 



j+s j+s+l 



= 0=A i+s A 



j+s n j+s+l ■ 



(2.15) 



Note that this is intended not as a statement on the op- 
erator algebra, but as a mere bookkeeping device for 
evaluating expectations. The assignment 

A] +s N j+s+1 =0 = N j+s A j+s+i (2.16) 
follows from Eq. ( 12.151 ). 

2. Right-exclusion map. The right-exclusion map de- 
scribed in Section UlAl is then implemented by making 
the substitution 



The assignment 

N j+s (t-Nj +s+ i) = n j+s 
follows from Eq. ( 12.171 ). 



(2.18) 



3. Re-indexing. Because the right-exclusion map in 
Eq. ( 12.17b merges the occupied site j + s and the empty 
site j + s + 1 to its right, operators to the right of site 
j+s+l must be re-indexed. The index j + s on the 
excluded chain becomes 



s'=0 



(2.19) 



A) +s (l-N j+s+1 )»a] +s 



(2.17) 



on the ordinary chain. Thus, two ending operators r 
sites apart in the /j-intervening-particle excluded expec- 
tation becomes r — p sites apart in the corresponding 
/^-intervening-particle ordinary expectation. 



III. CORRELATIONS IN THE BOSONIC AND 
FERMIONIC EXCLUDED CHAINS 

In this section, we make use of the tools developed in 
Sec. |n]to calculate three simple correlation functions within 
the ground states of the excluded chains of hardcore bosons 
and spinless fermions. In general, the intervening-particle 
expansion for excluded chain ground-state correlations must 
be evaluated numerically (even when each ordinary chain 
ground-state expectations in the sum can be expressed in 
closed form), keeping in mind a excluded chain with density 
N maps to an ordinary chain with density h — N 7(1 - N) (see 
Eq. Oh 

In Sec. [IITaI Sec. UlIBl and Sec. HiICl we show numeri- 
cal results for the two-point functions (B\Bj +r ) and (CjC/ +r ), 

and the four-point functions (NjNj+ r ), {B\B\,Bj +r By +r ) and 
{CjCj,Cj +r Cj' +r ) respectively. For the sake of easy ref- 
erence, we will call the two-point functions (BjBj +r ) and 

{C\Cj +T ) Fermi-liquid (FL) type correlation functions, even 
though their spatial structures depend on particle statis- 
tics. We will also call the four-point functions (NjNj +r ) = 
</>•;/»',/>•;.,/»',.,> , ^r;C ; C;. ,.(';.,•> CDW type correlation func- 
tions, and the the four-point functions (B\b\,B j +r B y +r ) and 

(C^C\,Cj +r C r+r ) SC type correlation functions. Both CDW 
and SC type correlations are identical for hardcore bosons and 
spinless fermions on the excluded chain, but the latter has the 
'superconducting' interpretation only for fermions. 

In Sec. IIII Al we will also explain in detail how nonlinear 
curve fits of the numerical correlations to reasonable asymp- 
totic forms as a function of r are done. Based on results from 
the nonlinear curve fits, we show how the Luttinger's theorem 
manifests itself, and how meaningful power-law exponents 
can be extracted. Similar analyses are done in Sec. IIII Bl and 
Sec. IIII CI as well as in Sec. [IV] for the three limiting ground 
states of the ladder model. 
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A. FL Correlations 



In the intervening-particle expansions of the two- 
point functions (BjBj +r ) and (CtC/ +r ), the nonvan- 
ishing terms map to /j-intervening-particle expecta- 
tions of the form {b\ \\ jp n jp ft jh (1 - n jh )b j+r - p ) and 

( c ) I\j p n j P rij,,(l _ n ih) c i+r-p)- B° tn can be evaluated in 
terms of two-point functions 



sinn7r|/ - j\ 
Ai ~ j\ 



(3.1) 



of the one-dimensional Fermi sea, after invoking the Jordan- 
Wigner transformation (see Appendix [A} for the former. 

As shown in the inset of Fig.[2] the FL correlation (B^Bj +r } 
was found to consists of a simple power law part, decaying 
with a smaller exponent ao, and an oscillatory power law part, 
decaying with a larger exponent a\. Multiplying {B\Bj +r ) by 
various simple powers of r, we find that yfr {BXB j +r ) asymp- 
totes to a constant with large r (as shown in the main plot of 
Fig. EJ, which suggests that ao = \- This is the correlation 
exponent predicted by Efetov and Larkin, in their study of the 
ordinary chain of hardcore bosons J£ 




FIG. 2: (Color online) Plot of ^fr (B\B j +r ) as a function of r for 
the particle densities N = 0.20 (red circles), N = 0.25 (green 
squares), N = 0.30 (blue diamonds). The colored curves shown 
are nonlinear curve fits of V? (B\B ; +r ) to the asymptotic form Ao + 

A[ r _ar i cos(2fc f r + (/>[). (Inset) Plot of (B^Bj +r ) as a function of r, 
showing that it consists of a simple power law part and an oscillatory 
power law part. 

Another important result comes from the unrestricted non- 
linear curve fits of ^{B'jBj +r ) to the asymptotic form Ao + 

A\ r ' cos(kr + (pi). In Fig. [3] we show a plot of the fitted 
wave number A; as a function of the density N of the excluded 
chain of hardcore bosons. As we can see, the fitted wave num- 
bers fall neatly onto the straight line k = 2kf = 2nN, where 
kf = nN is the Fermi wave number. The fact that kf ap- 
pears naturally in the numerical correlations is expected from 
Luttinger's theorem, which states that the volume of the re- 
ciprocal space bounded by the noninteracting Fermi surface 



is invariant quantity not affected by interactions, and applies 
in both Fermi and non-Fermi liquids i 2 ' 11,12 ' 13,14 - 15 ' 16 From this 
point onwards, we restrict the wave number of the oscillatory 
part of the correlation functions to kf, 2kf, or Akf in the non- 
linear curve fits. 




FIG. 3: Plot of the fitted wave number k (solid circles) against the 
density N of the excluded chain of hardcore bosons. The parameter 
k is obtained from the unrestricted nonlinear curve fit of -\fr <sjfi /+r ) 

to the asymptotic form A + A { r~ a i cos(kr + <j>{). The straight line is 
k F = 2nN for the Fermi wave number. 

From the restricted nonlinear curve fits, we find that A\ is 
large when a' { is large, and small when a' { is small. This sug- 
gests that neither of these parameters can be accurately deter- 
mined from our nonlinear curve fits, unless we further con- 
strain what values a[ can take. We also find that the quality 
of the nonlinear curve fit is good when is far from N - \, 
but deteriorates as we approach quarter filling. This suggests 
important physics in the FL correlation (B^Bj +r ) near quar- 
ter filling, which cannot be adequately accounted for by an 
asymptotic form Ao +A\ r~ a '> cos(kr + (pi). This loss of fit also 
affects the phase shift <f>\, presumably to a smaller extent, and 
the amplitude Ao of the simple power law, to an even smaller 
extent. These two parameters are plotted as functions of the 
density in Fig. [4] In the limit — > 0, we have essentially a 
dilute gas of hardcore (otherwise noninteracting) bosons, and 
thus (B^Bi+r) should include an overall factor of N. Thus we 



expect Ao 



as Af -> 0. In the limit N 



\, the excluded 



chain of hardcore bosons becomes increasingly jammed, and 
the relevant degrees of freedom are holes. For a dilute chain 
of holes, we expect {B\Bj +r ) to be proportional to the hole 

density, and thus Ao — » as N — » h. From our numerical 
results alone, it is hard to tell whether Ao reaches a maximum 
at = j (corresponding to a quarter-filled, n — |, ordinary 
chain of hardcore bosons) or = | (quarter-filled excluded 
chain of hardcore bosons). It is also hard to say anything defi- 
nite about the phase shift <p\, which might in fact be constant. 

In contrast, the FL correlation {C .Cj +r ) contains no simple 
power law part. A preliminary unrestricted nonlinear curve fit 
of this correlation to the asymptotic form A \ r~ a ' cos(nNr + 
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FIG. 4: Plot of the fitted amplitude Ao (top) of the simple power law 
part, and the fitted phase shift tp\ (bottom) of the oscillatory power 
law part, as functions of the density N of the excluded chain of hard- 
core bosons. 



4>\ ) suggests that a\ as 1 for all densities, i.e. may be universal 
a\ — 1 just as for noninteracting spinless fermions. How- 
ever, a more careful restricted nonlinear curve fit r {C\C ,- +r ) 
Aicos(nNr + <pi) show systematic deviations, as shown 
in Fig. [5] and therefore we perform an unrestricted fit to 
r{C-Cj +r ) = M r 1- " 1 cos(nNr + (pi). The fitted parameters 
A\, 1 - a\, and <f>\ are shown in Fig. [6] 

At very low densities N — > 0, our dilute chain of nearest- 
neighbor excluding spinless fermions will behave like nonin- 
teracting fermions, so we expect 



+- 




(C)C j+r ) 



sin nNr 



(3.2) 



FIG. 5: Plots of r(CtCj +r ) as functions of r (solid circles) for par- 
ticle densities N = 0.10 (top), N = 0.25 (middle), and N = 0.40 
(bottom). Nonlinear curve fits of r(C .Cj +r ) to the restricted asymp- 
totic form A\ cos(nNr + </>i) (solid curves) show systematic devia- 
tions, which can be accounted for by a unrestricted asymptotic form 
4j yfi-cri) eos(;riV> + fa). 



From our curve fits, we see that in this limit, A\ — > \jn — 
0.31831 .. ., 1 - a\ -> 0, and <f> x -> 3tt/2, and thus the FL 
correlation does indeed go to the low density limit in Eq. (13.2b . 

Also, in the half-filling limit N — > \, the chain become 
more and more congested, making it increasingly difficult to 
annihilate a spinless fermion at site j + r, find an empty site 
j to create a spinless fermion, without running afoul of the 
nearest-neighbor exclusion constraint. This tells us that A\ 
must vanish as — > i, which is hinted at in Fig. [6] However, 
the vanishing amplitude is only half of the story in this limit, 
the other half being the rate at which the correlation decay 
with increasing separation. In fact, very close to N = j, we 
expect the ground-state physics of the chain of rung-fermions 
with infinite nearest-neighbor repulsion to be describable in 
terms of a low density of holes. Naively, we would expect 
from such a low-density-of-holes argument that {C .Cj +r ) de- 
cay as r~'. Instead, the nonlinear curve fits at A^i < \ tells us 
fhatai < 1. 

Thinking about this nearly-half-filled limit more carefully, 
we realized that what we called 'holes' are really domain 
walls separating a region in which the spinless fermions sit 
on odd sites, from a region in which the spinless fermions sit 



on even sites. The FL correlation (C C, +r ), which can be writ- 
ten as a hole-hole correlation function, then depends on how 
many holes there are between j and j + r. The idea is that, 
in order to annihilate a hole (create a spinless fermion) at site 
j + r and create a hole (annihilate a spinless fermion) at site 
j, we must first find a configuration with a hole at j + r. Such 
a configuration will have spinless fermions at sites j + r - 2, 
j + r - 4, . . . , until we encounter another hole at j + r — Is, 
and then the sequence of spinless fermions will thereafter be 

at sites j + r — 2s - 1, j + r — 2s - 3, If r is even, (Cjc /+ ,-) 

receives nonzero contributions only from those configurations 
with an even number of intervening holes, whereas if r is odd, 
{C\Cj +r ) receives nonzero contributions only from those con- 
figurations with an odd number of intervening holes. This is 
very similar in flavor to the intervening-particle expansion of 
the two-point function (b^bj +r ) of a chain of ordinary hard- 
core bosons, except that (b\b; +r ) receives positive contribu- 
tions from configurations with an even number of intervening 
particles, and negative contributions from an odd number of 
intervening particles. Therefore, in the limit — > i, we find 
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FIG. 6: Plot of the fitted amplitude A, (top), the fitted exponent 
a i (middle), and the fitted phase shift 0i of the leading oscilla- 
tory power-law decay in the FL correlation {C\Cj +r ), as functions 
of the density N of the excluded chain of spinless fermions. In 
these plots, the solid curves are for fits to A] cos(7t/Vr + (p t ), 
whereas the dashed curves are for fits to A t r 1 ^" 1 cos(nNr + 4>i) + 
Ai r 1- " 2 cos[;r(l - N)r + <pi\. The dotted line in the bottom plot is a 
straight line from 5tt/4 at N = to 3tt/2 at N = j to guide the eye. 



B\ appears to be density-dependent, and is very close to being 
fr = 3 + f(s-#)- C 3 - 3 ) 




that the FL correlation (CjC /+ , ) maps to a string correlation 
of holes. Bosonization calculations then show that this string 
correlation of holes decay as a power law, with correlation 
exponent a i = i.— 



0.1 
0.0 
-0.1 



B. CDW Correlations 

Important physics can also be learnt from the nonlin- 
ear curve fitting of the CDW correlations (NjNj +r ) = 

</>•;/»',/>•,./>•;.,> = (C'lCiC'^Ci). First we tried to fit the sub- 

,/ j j +r j j j j'' j 

tracted CDW correlation to the asymptotic form (NjNj +r ) - 
(Nj){Nj +r ) = B\ r~ Pl cos(2nNr + 6\), but found the quality 
of fit deteriorates as N — > 0, as shown in Fig. [7] We under- 
stand this as follows: for N — > 0, the dimensionless quantity 
£ = Nr is small, and the poor fit indicates that {NjNj +r ) con- 
tains contributions from a term that decays more rapidly than 
B\ r~&' cos(27iNr + 9\). If we assume that this faster decay- 
ing term is a simple power law of the form Bo r ~^ 2 , and fit 
the CDW correlation to (NjN j+r ) = B x cos(2^7Vr + 9\) + 
Bo r~P 2 , we found the quality of the nonlinear curve fit is im- 
proved, after dropping data points r — 2, 3 from the fit, as 
shown in Fig. [7] We can include the simple power law decay 
term in the nonlinear curve fit throughout the entire range of 
^V, but the parameters Bo and B2 cannot be reliably determined 
beyond = j. Therefore, the parameters B\, f3\, and #i are 
the only parameters that can be reliably determined across the 
whole range of densities. From Fig.[8j we find that 6\ changes 
very little over the whole range of densities, and remains close 
to 7T/16. On the other hand, the leading correlation exponent 
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FIG. 7: Nonlinear curve fits of the subtracted CDW correlation 
(NjNj +r ) (solid circles) in the bosonic/fermionic ground states of the 
excluded chain, at densities N = 0.10 (top), = 0.25 (middle), and 
N = 0.45 (bottom). Above quarter filling (middle and bottoms plots), 
the subtracted CDW correlations can be fitted very well to the sim- 
ple asymptotic form B[ r^ 1 cos(2nNr + 6i) (solid curves), whereas 
at low densities (top plot), the subtracted CDW correlations deviate 
significantly from this simple asymptotic form. The nonlinear curve 
fit improves only after we add a simple power-law correction term 
Bi f^ 1 , giving the dashed curves. 

Just as for the FL correlation, we need to again think hard 
about the hole physics of the ground-state CDW correlation 
very close to half-filling. In this limit, the CDW correlation 
(NjNj +r ) can be written in terms of the hole-density-hole- 
density correlation {PjPj +r ), where Pj = 1 - Nj is the hole 
occupation number operator on site j. Using an argument sim- 
ilar to the one used for the FL correlation (Ctc,- +r ) near half- 

j J 

filling, we realized that if r is even, (PjPj +r ) receives nonzero 
contribution from configurations with an even number of in- 
tervening holes, whereas if r is odd, {PjPj +r ) receives nonzero 
contribution from configurations with an odd number of inter- 
vening holes. This means that our hole-density-hole-density 
correlation (PjPj +r ) must be mapped to a string correlation of 
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FIG. 8: Plot of the fitted amplitude Bi (top), the fitted exponent 
By (middle) and the fitted phase shift Q\ (bottom) of the leading 
oscillatory power law decay of the subtracted CDW correlation 
(NjNj +r ) - (Nj) (Nj+r), as functions of the density N of the excluded 
chain of hardcore bosons/spinless fermions. 



a chain of noninteracting spinless fermions. Bosonization cal- 
culations tell us that this string correlation decays as a power 
law with correlation exponent a = i consistent with the con- 
jectured behaviour, Eq. ( 13.31 ). at the special value N = 



C. SC Correlations 

In contrast to the FL and CDW correlations, the SC corre- 
lation (Aj__ 2 AjAj +r Aj +r+ 2) on the excluded chain has a rather 
more complex structure. The SC correlation is always nega- 
tive, and oscillations are highly suppressed, suggesting that it 
is the sum of a simple power law and an oscillatory power 
law. To improve the reliability of the nonlinear curve fits, 
we prescale the SC correlation by multiplying it by r 7 ^ 4 . 
This strange exponent is chosen because it is closest to the 
rate at which the simple power law decay for various densi- 
ties. After dropping data points for r = 2,3,4, good fits to 
the asymptotic form r 7/4 <A}_ 2 AU j+( .A j+ ,. +2 > = C r 7/4 ^° + 
Ci r 7/4 ^' cos(27T^Vr+^-i) were obtained. The nonlinear curve 
fits were improved marginally by adding a correction term of 
the form C 2 r 11 ^- cos(4nNr + X i) (see Fig. g). The fitted 
parameters are shown in Fig.[K)]as functions of the excluded 
chain density N. 

Unlike for the FL and CDW correlations, there are no ana- 
lytical bosonization calculations to help suggest values for the 
SC correlation exponents, so we used an ad-hoc process where 
we imposed trial values of the exponents, and let the nonlin- 
ear curve fitting program find the appropriate amplitudes and 
phase shifts. We found visually that the best fit of the numeri- 




-0.02 



FIG. 9: Nonlinear curve fits of the SC correlation 
r 7/4 (A._ 2 A^Aj +r Aj +r+ 2) (solid circles) in the bosonic/fermionic 
ground states of the excluded chain, at densities TV = 0.10 (top), 
N = 0.25 (middle), and N = 0.45 (bottom). After dropping data 
points from r = 2,3, 4, the SC correlations can be fitted well to the 



simple asymptotic form Co r 



.7/4-yo 



,7/4-n 



cos(2nNr + x\) (solid 



curves). The nonlinear curve fits are marginally improved by adding 
a correction term of the form C2 r 7/4 ~ y2 cos(4^/Vr + xi) (dashed 
curves). 



cal SC correlations appears to the mixed asymptotic form 

r 7/4 {A)_ 2 A)A j+r A j+r+2 ) = C' ri + 

-i - 
C\ r ? cos(2nNr + x\) + 

C 2 H cos(2nNr + X ' 2 ) + C' 3 ri. (3.4) 

IV. LADDER MODEL 

In this section we show how the analytical machinery devel- 
oped in Sec. [TT] can be adapted to calculate ground-state cor- 
relations in the ladder model of interacting spinless fermions 
given in Eq. ( 11.21 ). in three limiting cases where the ground 
states can be deduced from simple energetic arguments. An 
overview is given in Sec. I IV Al before we move on to de- 
tailed analyses and discussions of the three limiting cases in 
Sec. UVBl Sec. UVCl and Sec UVDl As with the chain mod- 
els, we assume that the ladder is finite, with j = 1 , . . . , L, and 
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FIG. 10: Plot of the fitted amplitudes Co and Ci (top), the fitted ex- 
ponents 70 and yi (middle), and the fitted phase shift x\ of the lead- 
ing simple power-law decay and the subleading oscillatory power- 
law decay in the SC correlation (A._ 2 A\Aj +r Aj +r +2)> as functions 
of the density N of the excluded chain of hardcore bosons/spinless 
fermions. 



subject each of its legs i — 1, 2 to open boundary conditions. 
Exact solution for the infinite ladder is then obtained by tak- 
ing L — > 00 keeping the particle density N2 fixed. Just as for 
the chain models, we expect in this limit that the ladder exact 
solutions would not depend on which boundary conditions we 
used. 



A. The Three Limiting Cases: An Overview 

For the ladder model described by Eq. ( 11.21 ). with V — > 00, 
the ground state is determined by the two independent model 
parameters, t±/t\\ and f /t\\, and the density i?2. For fixed N2, 
the two-dimensional region in the ground-state phase diagram 
is bounded by three limiting cases, 

(i) the paired limit t' » t\\, t±, which we will discuss in de- 
tail in Sec. IIVBI In this limit, we find SC correlations 
dominating at large distances (though, as for hardcore 
bosons, CDW correlations inevitably dominate at short 
distances). Based on our numerical studies in Sec. lIVBl 
the leading SC correlation exponent appears to be uni- 
versal, with a value of y — i while the leading CDW 
correlation exponent /3 is nonuniversal. In this limit, FL 
correlations are found to decay exponentially. A stag- 
gered form of long-range CDW order also appears; 

(ii) the two-leg limit t ± <k t\\, t' = 0, which we will dis- 
cuss in detail in Sec. IIVCI In this limit, the two legs of 
the ladder are coupled only by infinite nearest-neighbor 
repulsion. The dominant correlations at large distances 
are those of a power-law CDW, for which we find nu- 
merically to have what appears to be an universal cor- 
relation exponent of /? = i. In this limit, the leading 
SC correlation exponent was predicted analytically to 



be y = 2, while FL correlations are found to decay ex- 
ponentially; 

(iii) the rung-fermion limit t± » t\\, t' = 0, which we 
will discuss in detail in Sec. IIVDI In this limit, the 
particles are effectively localized onto the rungs of the 
ladder. When the ladder is quarter-filled, a true long- 
range CDW emerges in the two-fold degenerate ground 
state. Below quarter-filling, we find numerically that 
the CDW power-law correlation dominate at large dis- 
tances, with a leading non-universal correlation expo- 
nent p — 5 + § (5 _ ^Vi)- The leading FL correlation ex- 
ponent was also found numerically to be non-universal, 
with values going from a — 1 to a — |. The SC corre- 
lation exponent, on the other hand, was found numeri- 
cally to be universal, with value y =\- 

To zeroth order (i.e. without plunging into first-order per- 
turbation theory calculations), the ground-state phase diagram 
can be obtained by interpolating between these three limit- 
ing cases. There will be three lines of quantum phase tran- 
sitions or cross-overs, which at quarter-filling, separate the 
long-range CDW (LR-CDW), power-law CDW (PL-CDW), 
and SC phases. At quarter-filling, we expect these three lines 
of critical points or cross-overs to meet at a point on the phase 
diagram. If we have three lines of true critical points, this 
point would be a quantum tricritical point. We therefore end 
up with a ground-state phase diagram which looks like that 
shown in Figure QTJ 



t±/t\\ 



LR-CDW 



SC 



1 

PL-CDW 1 
1 
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FIG. 11: The zeroth-order ground-state phase diagram of the ladder 
model given by Eq. dl.2t . The three limiting cases we can solve 
exactly are shown as the two dots (cases (ii), power-law CDW (PL- 
CDW) and (iii), long-range CDW (LR-CDW)), and the thick solid 
line (case (i), SC). 



B. The Paired Limit 

In this subsection, we solve for the ground-state wave func- 
tion, and calculate various ground-state correlations in the 
paired limit t' » t\\,t ± . In this limit, the Hamiltonian in 
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Eq. ( 11.2b simplifies to 

H fV = -t' X X (<■;.,''/• I./- 1 <v.;o + '';,.»m.mo,) 

-''XX (';•:./''•/• | ; - + c I + ij+2 n W+i c i+i.j) 

+ v X X Vy+i + y X X i.j- 
> j > j 

(4.1) 

In Sec. II V B 11 we explain how pairs of spinless fermions 
are bound by correlated hops in this limit, and the degrees 
of freedom in the system become mobile bound pairs with 
infinite nearest-neighbor repulsion. These bound pairs come 
in two flavors, determined by the specific arrangement of the 
two bound-pair particles around a plaquette. These flavors 
are conserved by correlated hops if the length of the ladder is 
even, and hence the ladder ground state is two-fold degener- 
ate. We then describe how these two degenerate ladder ground 
states can be mapped to a excluded chain of hardcore bosons, 
then to an ordinary chain of hardcore bosons, and finally to a 
chain of noninteracting spinless fermions. 

In Sec. IIV B 21 we calculate the SC and CDW correla- 
tions, using the intervening-particle expansion described in 
Sec. Ill B I We then use a restricted-probability argument in 
Sec. irVB~3l to show that FL correlations decay exponentially 
with distance, governed by a density-dependent correlation 
length, in this paired limit. We find, as expected from mak- 
ing the absolute correlated hopping amplitude t' large, that 
SC correlations dominate at large distances. 



1. Bound Pairs and Ground States 

In the paired limit f » £||,/_l, we solve for the ground 
state of the simplified Hamiltonian given by Eq. d4.U , which 
admits only correlated hops. Because of this, isolated spin- 
less fermions cannot hop at all; by contrast, a pair occupy- 
ing diagonal corners on a plaquette can perform correlated 
hops. Therefore, for an even number of spinless fermions, 
ground-state configurations consist of well-defined bound 
pairs, which are effectively bosons. We say that a bound pair 
at (1,7) and (2, j+ 1) has even (resp. odd) flavor if its two sites 
are even (resp. odd) sites. In this limit of f /t\\,t'/t± — > oo, a 
particle on rung j can only hop to rung j ± 2. This moves the 
bound pair's center of mass by one lattice constant, without 
changing its flavor. The degrees of freedom in this limiting 
case thus becomes bound pairs with definite flavors hopping 
along a one-dimensional chain. 26 We write these hardcore bo- 
son operators in terms of the spinless fermion operators as 

+ \c\ f cl j odd; 
B\, - I \' J 2 ' } l v J (4.2) 

[ C ij+Aj> J even > 

and 

t \c\ ;j -i c i ;> j °dd; 
B\_ = { b J+ \ 2 ' J (4.3) 

;, ~ K/L+i' J even - 



where we order first with respect to the leg index, and then 
with respect to the rung index of of the ladder. 

Since bound pairs cannot move past each other along the 
chain, the P-bound-pair Hilbert space breaks up into many 
independent sectors, each with a fixed sequence of flavors. 
The P-bound-pair problem in one sector is therefore an in- 
dependent problem from that of another P-bound-pair sector. 
The minimum energy in each sector can be very crudely de- 
termined by treating the P-bound-pair problem as a particle- 
in-a-box problem, where each bound pair is free to hop within 
a 'box' demarcated by its flanking bound pairs. 
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FIG. 12: The closest approach two bound pairs can make to each 
other, if (a) they both have even flavors; (b) they have opposite fla- 
vors; and (c) they both have odd flavors. 

As shown in Figure [12] two bound pairs with the same fla- 
vor can get within a separation r — 2 of each other, whereas 
two bound pairs with different flavors can only achieve a clos- 
est approach with separation r — 3. Therefore, for a fixed sep- 
aration between flanking bound pairs, the kinetic energy of the 
'boxed' bound pair is lowest when all three bound pairs have 
the same flavor. Repeating this argument for all bound pairs, 
we realized therefore that the two-fold degenerate ground state 
lies within the all-even and all-odd sectors. In these sectors, 
bound pairs cannot occupy nearest-neighbor plaquettes, i.e. 
we are dealing with an excluded chain of hardcore bosons^ 

The twofold degeneracy between all-even and all-odd sec- 
tors represents a symmetry breaking with long-range order of 
a staggered CDW type (in terms of fermion densities). It 
may be viewed as breaking the invariance under reflection 
about the ladder axis of the original Hamiltonian as given in 
Eq. ( 11.2b . 28 Thus the quantum-mechanical problem of a ladder 
with density N2 is mapped to the quantum-mechanical prob- 
lem of an excluded chain with density N - N2- 



2. SC and CDW Correlations 

Three simple correlation functions, the FL, CDW, and SC 
correlations, were computed for the excluded chain of hard- 
core bosons in Sec. [TTlJ On the ladder model in the paired 
limit, these correlations must be interpreted differently. In 
mapping the ladder model to the excluded chain, we re- 
place a pair of spinless fermion by a hardcore boson opera- 
tor, i.e. c\ .ct .,, — > B\. Thus the FL correlation (B\b ,•+,) of 

lj 2j+l j j J +r > 
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hardcore bosons actually corresponds to a SC correlation of 
the fermion model on the ladder. Depending on which of the 
two degenerate ladder ground states we are looking at, the SC 
operators are 



V2 
1 

V2 



t t , t t \ 

■' 1 T c 1,7+1 L 2,y' 



2J+1 



At = -i=(_iy( c t 



C 2,;+l C lj+l l -2J 



Cl :), 



(4.4) 



such that 



{A^A jw > s = <A;A jt „>.. 



0. 



(4.5) 



Because of Eq. d4.5l ), we shall drop the indices § and u from 
here on. From Sec. [Till we know that (AjA J+ ,} decays with 
separation r asymptotically as the sum of a simple (leading) 
power law and an 2A:/r-oscillatory (subleading) power law. 
The leading correlation exponent has been determined to be 
jo - \, while the subleading correlation exponent j\ cannot 
be reliably determined. 

The simplest CDW correlations are 



< C l,j C hj C l J+r C hj+r) , {c\fXjcl jJr C % j +T ) , 



(4.6) 



which we call the CDW-cr correlations. These are not easy to 
calculate, because they cannot be written simply in terms of 
the expectations of hardcore boson operators. In contrast, the 
CDW-7T correlations^ 



{B\B jB\ +r B j +r ) = {NjN j+r ) 



(4.7) 



can be evaluated with the help of the intervening-particle ex- 
pansion in Eq. ( 12.141 ). This was done in Sec. [Hi] where 
we found the subtracted CDW-;r correlation {NjNj +r ) - 
(Nj) (Nj+r) decaying asymptotically with separation r as a 
simple power law Bi r^ 1 cos(2£pr + 6\), with a nonuniver- 
sal leading correlation exponent [3\ — | + | (| - N2), and a 
universal phase shift of 6\ — 7r/16. 



3. FL Correlation: Explanation of Exponential Decay 

Unlike the SC and CDW-7T correlations, the FL correlations 
cannot be calculated easily in this paired limit, because the op- 
erators involved cannot be written in terms of hardcore boson 
operators. Nevertheless, we can still calculate it by making 
use of the fact that this correlation is very close to being the 
probability of finding a restricted class of configurations in the 
ground state. We then make use of the scaling form reported 
in Reffl8l to calculate the probability analytically. This idea is 
exploited again in Sec. lIV C~2l 

In this paired limit, the ground state consists exclusively 
of a superposition of bound pair configurations. Therefore, if 
we annihilate a spinless fermion on leg i, we must create an- 
other on the same leg elsewhere, and thus the only nonzero 

FL correlations are of the form (c! c, ,- +r ). In fact, to start with 

i >j ,J 



a paired configuration and end up with another paired con- 
figuration, after annihilating a spinless fermion at j + r and 
creating a spinless fermion at j, the initial and final configura- 
tions must contain a compact cluster of pairs between rung j 
and rung j + r, as shown in Fig. [13] 
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FIG. 13: Compact p-bound-pair cluster configurations making 
nonzero contribution towards the expectation of the FL operator 
product c[ .Cij + 2 P , which annihilates a spinless fermion on the right 
end of the compact p-bound-pair cluster, and creates a spinless 
fermion on the left end of the compact p-bound-pair cluster. "P, and 
*?/ are the ground-state amplitudes of the initial and final configura- 
tions respectively. 

Based on this compact cluster argument, we know that 
(cJjCij+r) = when r is odd. When r — 2p is even, 



(C}fi, j+ 2 P ) = Yj^i 



(4.8) 



<;,/) 



receives contributions from all pairs of configurations with a 
compact p-bound-pair cluster between the rungs j and j + 2p. 
Clearly, these products of amplitudes will depend on where 
the other bound pairs are on the ladder. However, if the ladder 
is not too close to half-filling, we expect *¥/ ~ *P,, so that on 
an infinite ladder, (cJjCij+2 P ) is very nearly the probability of 
finding a compact /j-bound-pair cluster, 10 



A^2 

{NjNj+2 ■ ■ ■ Nj+2p) = — (n jti j+i---n j+p ) 

= -idetG c (p), 
n 



(4.9) 



after using the relation ( 12. lib between excluded and ordi- 
nary expectations, where n = #2/(1 - N2) is the density 
of the ordinary chain. Here det Gc(p) is the determinant 
of the noninteracting-spinless-fermion cluster Green-function 
matrix Gc(p) for a cluster of p sites, which we can write as 1! ^ 



detGc(p) = ]~[* = r 



1 



+ r 



(4.10) 



where Ai are the eigenvalues of the cluster Green-function ma- 
trix Gc(p), and <pi are the single-particle pseudo-energies of 
the cluster density matrix pc, for the cluster of p sites in an 
infinite chain of noninteracting spinless fermions. 

For p » 1, Gc(p) has approximately hp eigenvalues which 
are almost one, and approximately (1 —n)p eigenvalues which 
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are almost zero, The determinant of Gc(p) is thus determined 
predominantly by the approximately (1 - h)p eigenvalues 
which are almost zero. For these Ai, e Vl » 1, and thus 



' i F +a-n)p \ 

fetG C {p)*Y\ e ^' = ex p - Z W 



(4.11) 



where Z/r is such that tpy = 0. Converting the sum into an 
integral, and using the approximate scaling formula in Reffl8l 
we find that 



det G c (p) ~ exp -p 



Jo 



f(n, x) dx 



(4.12) 



i.e. the probability of finding a compact p-bound-pair cluster 
decays exponentially with p in the limit of p » 1 . 

With this simple compact cluster argument, we conclude 
that the ladder FL correlation (cl jCij+ r ) decays exponentially 
with separation r as 



(cJ.c u+r >~exp[-r/f(iV 2 )], 
with a density-dependent correlation length 

2 



t(N 2 ) = 



f 

Jo 



l-n(2V 2 ) 



(4.13) 



(4.14) 



f(h(N 2 ),x)dx 



in the strong correlated hopping limit. From ReflTl we know 
that the scaling function f(h, x) depends only very weakly on 
n, and thus, at very low ladder densities N2 — > 0, the correla- 
tion length ^(A^) attains its minimum value of 



J f(0,x)dx 



(4.15) 



and the FL correlation (c. Cy +r ) decays most rapidly in this 
regime of N2 — > 0. This is expected physically, since a long 
cluster of occupied sites is very unlikely to occur at very low 
densities, with or without quantum correlations. 

In the regime of N2 — > 5, we find h — > 1, and thus the cor- 
relation length f (N) diverges according to Eq. (14. 14b . This di- 
verging correlation length tells us nothing about the amplitude 
of the FL correlation. Indeed, when the ladder becomes half- 
filled, the two degenerate ground states are inert bound-pair 
solids. Each of the half-filled-ladder ground-state wave func- 
tions consists of a single configuration whereby all available 
plaquettes are occupied by a bound pair, and it is not possible 
to annihilate a spinless fermion at the (j + r)th rung and create 
another at the jth rung. The FL correlation (cj jCi,j+ r ) is thus 
strictly zero in this half-filled-ladder limit. 



C. The Two-Leg Limit 

This subsection concerns the ground state in the two-leg 
limit t ± <K fy, t' = 0. Based on energetic considerations, 
we argue in Section IIVC II that there will be two degenerate 



ground states, within which successive spinless fermions are 
on alternate legs of the ladder. We call these the staggered 
ground states, and write their wave functions in terms of the 
Fermi sea ground-state wave function with the help of a stag- 
gered map between ladder configurations and ordinary chain 
configurations. We then calculate various ground-state corre- 
lations in Sec. IIVC 21 Sec. IIVC 31 and Sec. IIVC 41 where we 
show that the non-vanishing FL correlations decay exponen- 
tially with distance, governed by a density-dependent correla- 
tion length, while the CDW and SC correlations decay with 
distance as power laws. We find in this two-leg limit that the 
antisymmetric CDW correlation dominates at large distances. 



1. Ground States 

In the limit of fj_ — » 0, each spinless fermion on the two- 
legged ladder carries a permanent leg index, and thus the num- 
ber of spinless fermions P, on leg i are good quantum num- 
bers. Furthermore, successive spinless fermions along the lad- 
der cannot move past each other, even if they are on different 
legs, because of the infinite nearest-neighbor repulsion act- 
ing across the rungs. Consequently, the Hilbert space of the 
P-spinless-fermion problem breaks up into many independent 
sectors, each with a fixed sequence of leg indices. The P- 
spinless-fermion problem in one such sector is therefore an 
independent problem from that of another P-spinless-fermion 
sector. Noting that the closest approach between two particles 
on the same leg is r — 2, whereas that between two particles on 
different legs is r — 1, we invoke the same "particle-in-a-box" 
argument used for the paired limit in Sec. IIV B II to find the 
ground state for P spinless fermion on a ladder of even length 
L to be in a staggered sector, where successive particles are 
on different legs. There are two such sectors in a ladder with 
open boundary conditions, which we call sector 1 when the 
first fermion (from the left) is on leg 1, or sector 2 when it is 
on leg 2. 

Evidently this is a twofold symmetry breaking. (The bro- 
ken symmetry is that of reflecting the configuration about the 
ladder axis, which is a valid symmetry within the staggered 
sector.) This state has a form of long range order, in that the 
flavor alternates; however, that cannot be represented by any 
local order parameter, but only by a "string" order parame- 
ter^i 

Let us write |Vfi) and ^2) for the ground states in sectors 
1 and 2, respectively. A staggered configuration of P ladder 
spinless fermions in sector 1 can be mapped to a chain of P 
noninteracting spinless fermions using the staggered map 



. c; .. •• - c; ._ c . |0)i adder 



'lj>_i 2 ,j P 1 



C ■ C . 

Jl J2 



•eJ M C/,|0>ctata. (4-16) 



Using the same formula, with an exchange of leg index 1 <-> 2, 
configurations in sector 2 are similarly mapped^ 

Because the staggered map maps a ladder with density ^ 
onto a ordinary chain with density n - 2N2, we want a lad- 
der observable <9i a ddei an d its corresponding chain observable 
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Ochain to be such that 

(Oladder)ladder = \ (0chain)chain • (4.17) 

This is analogous to Eq. ( 12.1 11 ), which we derived when we 
map from a excluded chain to an ordinary chain. We use the 
subscripts 'ladder' and 'chain' just this once to distinguish be- 
tween ladder and chain expectations. This notation is cumber- 
some, so we will not use it again. Whether an expectation is a 
ladder expectation or a chain expectation should be clear from 
the context. 



2. FL Correlations: Exponential Decay 

Having solved the staggered ground states in terms of the 
one-dimensional Fermi sea, we calculate the FL, CDW, and 
SC correlations. There are four FL correlations at range r, 

(c[ .c j+r ), (c\jC2,j+ r ), <4j c U+r> and (c ] 2 j c 2 j +r ). From the 
staggered nature of \ y ¥±), we know that 

(c\ jClJ+r) = <Cj f2,j+r) \ 

: :' (4.18) 

(c[jC 2 J + r) = {C\f\,j + r) 

in both ground states. The inter-leg FL correlations vanish, 
i.e. 

(c\ f%j + r) = = (cljCij+r) , (4.19) 

because annihilating a particle on one leg and creating a par- 
ticle on the other leg disrupts the stagger configuration. 
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FIG. 14: Annihilation of a spinless fermion at site (I, j + r), followed 
by creation of a spinless fermion at site within a staggered 

ground-state configuration leads to a staggered ground-state config- 
uration, when there are no intervening particles between rungs j and 
j+r. 

The intra-leg FL correlation (c^c 1J+ ,), which is nonzero, 
receives contributions only from initial and final staggered 
configurations in which there are no intervening particles be- 
tween rungs j and j + r, for example, those shown in Fig. [14] 
This tells us that 

<cljC U+r ) = \ (ct(l - n j+l ) •■■(!- « /+ ,-i)c >+ ,.> , (4.20) 

when we map the ladder model to the chain model. This cor- 
relation is evaluated numerically, and shown in Fig.[l5j where 
we see the staggered ground-state FL correlations decaying 



exponentially with separation r. This asymptotic behaviour 
can again be understood using a constrained probabilities ar- 
gument similar to that used in Sec. |IVB2l except that instead 
of a compact cluster, the relevant probability P{r) is that of 
finding a gap at least r in length within the one-dimensional 
Fermi-sea ground state. 
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FIG. 15: The infinite-ladder FL correlations (c. .Cij +r ), i = 1,2, as a 
function of the separation 1 < r < 15 for ladder densities N 2 = 0.20, 
0.25 and 0.30, in the two-leg limit tjt% -* 0, f = 0. 

Applying a restricted probability argument similar to the 
one outlined in Sec. lIVB31 we know this probability is simply 
the zero-particle weight 

P{r) = w = det(l - Gd>)), (4.21) 

of the density matrix of a cluster of r contiguous sites in the 
chain of noninteracting spinless fermions. For r » 1, the clus- 
ter Green-function matrix Gc(r) has approximately (1 - n)r 
eigenvalues which are almost zero, and hr eigenvalues which 
are almost one. The determinant of 1 - Gc(r) is thus essen- 
tially determined by the approximately hr eigenvalues which 
are almost one. Using this fact, we calculate the asymptotic 
form of P(s) to be 

P(r) « exp j-r J /(l - n, x) dxj , (4.22) 

wh ere f (», x) is the universaling scaling function identified in 
Refill Eq. d4~22b expl ains the observed exponential decay 
of (ctjCjj+r) in Fig. [15] We note further that as « — > 1 (or 

equivalently, N 2 — > 5), the FL correlations decay fastest ex- 
ponentially, whereas as h — > (equivalent to N 2 0), the 
exponential decay is the slowest. We expect these behaviours 
physically, because it is more likely to find a long empty clus- 
ter when the density is low, and less likely to find a long empty 
cluster when the ladder is closed to half-filled. 
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3. CDW Correlations 

Next, we calculate the CDW correlations, for which the 
four simplest at separation r are, 



( C \,j C lJ C \,j + r C lJ+r} = {c[fXf' 2 j +r C 2 j +r ) , 
( C lj C lJ C 2,j+r C 2J+r) = (4,j C 2>j C h+r C l.J+r) ■ 



(4.23) 



Because of the staggered nature of the ground states, config- 
urations making nonzero contributions to («;,,«,■,/+,) are those 
which map to noninteracting spinless fermion configurations 
in which the sites j and j + r are occupied, with an odd number 
of intervening particles between them. Similarly, configura- 
tions making nonzero contributions to (n,jn,vj +r ) , i + i', are 
those which map to noninteracting spinless fermions in which 
the sites j and j + r are occupied, with an even number of 
intervening particles between them. 
Defining the density operators 



n±j = n Uj ± nzj 



(4.24) 



which are symmetric and antisymmetric with respect to reflect 
along the ladder axis, we find that 



and 



(n + jn-j +r ) = = (n-jn + j +r ) , 



(n + jn + j +r ) = (njUj+r) = E+(r), 



(4.25) 



(4.26) 



which we call the CDW+ correlation. This is identical to the 
CDW correlation of the one-dimensional Fermi sea, which we 
know decays as an oscillatory power law 

(n+jn + j +r ) - (n+j) {n + j +r ) ~ r 2 cos(2k F r). (4.27) 

There is also the CDW— correlation, 

(n-jn-j+ r ) = 2 (<«i. )«!,/+,) - (nijn 2 j +r yj = £-(r), (4.28) 

associated with n-j. This is identical to the subtracted CDW- 
correlation, since («-/) = («ij-«2j) = 0. Evaluating this 
expectation numerically, we find that at all ladder densities 
N2, S-(i?2, r) oscillates about a zero average with wave vector 
2kf , and a decaying amplitude. A preliminary unrestricted 
nonlinear curve fitting to the asymptotic form 



S_(AT 2 , r) = B r~ p o + B l r~ P] cos(2k F r + 6y), 



(4.29) 



where B[ r~P' cos(2kfr + 6{) is the leading asymptotic be- 
haviour, and Bo is a correction term, suggests that the 
leading correlation exponent may actually be universal, taking 
on the value B\ — j. Further nonlinear curve fitting, restrict- 
ing Bi = |, tells us that only the parameters B\ and 9\ of the 
leading asymptotic term can be reliably determined. These 
are shown in Fig. [16] From this restricted curve fit, it appears 
that the phase shift might also be universal, taking on value 
01 = n. 
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FIG. 16: Plot of the fitted amplitude B[ (top) and fitted phase shift 
6i (bottom) of the leading oscillatory power-law decay, as functions 
of the ladder density N2, for the CDW- correlation in the staggered 
ground state of the ladder model, in the two-leg limit t±/t^ — > 0, 
i = 0. 



4. SC Correlations 
The simplest SC correlations at separation r are 

( C l,/2j+l C W+'' C 2,;+ r +l> = (c[ jC\ j+l C 2 J +r Cij +r+ i} , 



( C l,j C 2J+\ C2 J+r C hj+r+l) ~ { C 2J c \J+\ c l,j+rC2J+r+l) ■ 



(4.30) 



Correlations of the type (c]jcj, j +l Cij+ r Ci>j+ r +i) receive 
nonzero contributions from configurations containing an even 
number of intervening particles between rungs j + 1 and j + r, 
whereas correlations of the type (c! ; .c!, . +1 c/<j +/ -c,-j +r+ i) re- 
ceive nonzero contributions from configurations containing an 
odd number of intervening particles between rungs j + 1 and 
j + r. Defining the paired operators 

which are symmetric and antisymmetric with respect to reflec- 
tion about the ladder axis, we find that 



(A; j .A_ J+r > = = (Al J A +J+r >, 
and that the SC+ correlation 

(Al A +J+r ) = <c]c] +1 c j+r c j+ r + i) = n+(r) - 



at large separations. 
The SC- correlation 



<Al J A_ J+r > = II_(r) 



(4.32) 



(4.33) 



(4.34) 



must be evaluated numerically. We find that, just like 2_(r), 
Tl-(r) oscillates about a zero average with wave vector 2kf, 
and a rapidly decaying amplitude. To improve the quality of 
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the nonlinear curve fitting, we fit r 2 II_(r) to the asymptotic 
form 

r 2 EL(r) = C r 2 ^ + Ci r 2 ^' cos(2k F r + X i), (4.35) 

where C\ r 2 ^' cos(2kfr + x\) is the leading asymptotic be- 
havior, while Co r 2 ~P° is a correction term. A preliminary un- 
restricted fit suggests that the leading correlation exponent is 
universal, and takes on value ySi = |, Further restricted non- 
linear curve fitting tells us that only the parameters C\ and;^ 
can be reliably determined. These are shown in Fig. [17] where 
we see that the amplitude C\ exhibits symmetry about quar- 
ter filling, which is a kind of particle-hole symmetry, and that 
the phase shift x\ — 71 [l + \ {j ~ ^2)] is non-universal, but 
depends linearly on the density $2. 
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FIG. 17: Plot of the fitted amplitude Ci (top) and fitted phase shifty 
(bottom) of the leading oscillatory power-law decay, as functions of 
the ladder density N 2 , for the SC- correlation in the staggered ground 
state of the ladder model, in the two-leg limit t x /t\\ — > 0, f = 0. 



D. The Rung-Fermion Limit 

In this subsection, we look at the rung-fermion limit t ± » 
f||, t' = 0. We argue in Section ITVD II that in this limit, each 
spinless fermion spends most of its time hopping back and 
forth along the rung it is on, and only very rarely hops along 
the legs to an adjacent rung. Therefore, each spinless fermion 
will be in a quantum state very close to the symmetric eigen- 
state of one rung, and we can think of the ladder of spinless 
fermions with density N2 in this limit as essentially an ex- 
cluded chain of fermions with density N = 2Nz- For N2 < \, 
the ground state of this excluded chain of rung-fermions has 
been solved in Sec. Ill Al The FL, CDW, and SC correlations 
have also been calculated in Sec. [Till so we will not repeat 
them here. 

At N2 = t, the ground state is a 'dynamic solid' phase, 
in which rung-fermions occupy either all the even rungs, or 
all the odd rungs, and cannot hop along the legs to adjacent 



rungs because of the infinite nearest-neighbor repulsion be- 
tween them. For N2 > |, we describe in Sec. IIV D~2l how 
the system will phase separate into a high-density inert solid 
phase, in which spinless fermions cannot hop at all, and the 
lower-density 'dynamic solid' phase. In this phase separation 
regime, the FL, CDW, and SC correlations cannot be calcu- 
lated. 



1. Ground States 

In the limit of t\\/t x — > 0, a spinless fermion spends most of 
its time hopping back and forth along a rung, and only very 
rarely hops along the leg to an adjacent rung, where it will 
spend a lot of time hopping back and forth, before hopping 
along the leg again. Because of this long dwell time on a 
rung, the spinless fermion is very nearly in the rung ground 
state 



c>>. 



(4.36) 



Let us call a spinless fermion in the rung ground state a rung 
fermion in short. Rung-fermions inherit the infinite nearest- 
neighbor repulsion of the bare spinless fermions, and there- 
fore two rung-fermions in adjacent rungs experience infinite 
nearest-neighbor repulsion as well. With this insight, we find 
that the full many-body problem of spinless fermions with 
infinite nearest-neighbor repulsion on the two-legged ladder 
with density N2 reduces to the problem of an excluded chain 
with density N = 2N2 of spinless rung-fermions. 

The latter problem was solved in Sec.[n]and Sec.[in]for ex- 
cluded chain densities N < i. In the special case of quarter- 
filling on the ladder, N2 — \, spinless fermions occupy alter- 
nate rungs. These are free to hop along the rungs that they 
reside on, but cannot hop along the legs, for non-vanishing 
values of t\\/t±. Even virtual processes in which a spinless 
fermion on rung j hops along the leg to an adjacent rung and 
back are essentially forbidden by the infinite nearest-neighbor 
repulsion, because such virtual processes, which has a time 
scale of 0(l/f||), would not be complete when the spinless 
fermion on the next-nearest-neighbor rung hops across the 
rung, which occurs on a time scale of 0(1 /t±). Virtual pro- 
cesses such as these only become energetically feasible when 
the two time scales become comparable, i.e. when t\\ < fj_. 
Therefore, over a wide range of anisotropics t\\/t x , the spinless 
fermions in the quarter-filled ladder with tf — can hop back 
and forth along the rungs they are on, but cannot hop to the 
neighboring rungs. This gives rise to a symmetry breaking, 
where the spinless fermions are either all on the even rungs, 
or they are all on the odd rungs. Because translational symme- 
try along the ladder axis is broken in the quarter-filled ladder 
ground states, we think of these as 'dynamic solids', since the 
constituent spinless fermions are constantly hopping back and 
forth along the rungs. In this limit, the only non-vanishing 
correlation is the rung-fermion CDW correlation 



^ , r even; 
(NjN J+r )= 2 ' 

0, rodd, 



(4.37) 



16 



i.e. there is true long-range order in the quarter-filled ladder 
ground state in the limit of t± » fy, t' — 0. 



2. Phase Separation 

In this rung-fermion limit t ± » t\\, t' = 0, the system phase 
separates for ladder densities N2 > 4. As shown in Fig. [T8l 
when the ladder is above quarter-filling, some of the spinless 
fermions will go into a high-density inert solid phase with 
density N2 - where spinless fermions are arranged in a 
staggered array, and therefore cannot hop at all. These spin- 
less fermions contribute nothing to the ground-state energy. If 
fy is comparable to t±, the rest of the spinless fermions will go 
into a fluid phase, whose density is N2 < 4. These spinless 
fermions are free to hop back and forth on the rungs they are 
on, and occasionally to the neighboring rungs, when permit- 
ted by nearest-neighbor exclusion. These contribute a density- 
dependent total kinetic energy to the ground-state energy. The 
ground-state composition depends on whether the kinetic en- 
ergy gained per particle, by removing a spinless fermion from 
the solid phase and adding it to the fluid phase, outweighs the 
decrease in kinetic energy per particle that results from the 
fluid becoming more congested. 
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FIG. 18: Phase separation of a greater-than-quarter-filled ladder 
of spinless fermions with infinite nearest-neighbor repulsion into 
a high-density inert solid phase (immobile spinless fermions) with 
N2 = 5, and a low-density fluid phase (mobile spinless fermions 
shown with arrows) with N 2 = \, in the rung-fermion limit t± » f||, 
f = 0. 

When t± becomes large compared to fy, which is the limit 
we are interested in, it becomes energetically favorable, al- 
ways, to remove one spinless fermion from the inert solid 
phase, and add it to the fluid phase, if its density is N2 < 7. 
This is because the kinetic energy penalty to make the fluid 
becoming more congested, which is of <9(f||), is more than 
compensated for by the kinetic energy gain of t ± for an ex- 
tra spinless fermion freed to hop back and forth along a rung. 
Iterating this argument, we find then that, for t ± » fy, and the 
overall density N%> h, the system will phase separate into an 
inert solid phase with density N2 — h, and a dynamic solid 
phase with density N2 = j- For example, if the overall den- 
sity is M = I > 7, we will find that i of the total number 
of spinless fermions will be in the inert solid phase, while the 
other I of the total number of spinless fermions will be in the 
dynamic solid phase. 



V. SUMMARY AND DISCUSSIONS 

In this paper, we established a one-to-one correspondence 
between f-particle configurations on the excluded chain and 



P-particle configurations on the ordinary chain using the right- 
exclusion map. We then showed that the Hamiltonian matrices 
of the models given in Eq. ( 11.11 ) and Eq. ( 11. 3t are identical, 
therefore solving for the ground states of the former in terms 
of those of the latter. These results were obtained for finite 
chains subject to open boundary conditions, but continues to 
hold for infinite chains. 

Based on this one-to-one correspondence between ground 
states, we showed that the ground-state expectation {O) of an 
excluded chain observable O can be evaluated using Eq. ( 12.1 11 ) 
in terms of the ground-state expectation (O') of a carefully 
chosen corresponding observable O' on the ordinary chain. 
We then developed the method of intervening-particle expan- 
sion, to write the ground-state expectation {OjOj +r ) of a prod- 
uct of local excluded chain operators Oj and Oj +n first as 
a sum over excluded chain expectations (OjO p Oj +r ) condi- 
tioned on the occupations of the sites between j and j + r, and 
then as a sum over the corresponding ordinary chain expecta- 
tions (0'jO' p O' j+r _ p ). 

Using these analytical results from Sec. HH we calculated 
the FL, CDW, and SC correlations of the excluded chains of 
hardcore bosons and spinless fermions in Sec. [Till Based on 
nonlinear curve fits of the numerically evaluated correlations, 
to reasonable asymptotic forms, we find all three types of cor- 
relations decaying with separation r as power laws, for hard- 
core bosons as well as for spinless fermions. More interest- 
ingly, we find for both hardcore bosons and spinless fermions 
a universal exponent J\ - \ for the oscillatory power-law 
decay of the SC correlation, but a non-universal, density- 
dependent, exponent ySj = \ + 5(5 — N) for the oscillatory 
power-law decay of the CDW correlation. Also, the leading 
asymptotic behaviour for the hardcore boson FL correlation 
was found to a non-oscillating power-law decay with univer- 
sal exponent ao = \, while that for the spinless fermion FL 
correlation was found to be oscillations in a power-law enve- 
lope, with a non-universal exponent that approaches a\ = 1 as 
N -> 0, andffi = ± as N -» \. 

We then analyzed our spinless-fermion ladder model, 
Eq. ( 11.21 ). in Sec. [IV] This ladder model can be solved ex- 
actly in three limiting cases: (i) the paired limit t' » fy, t ± ; (ii) 
the two-leg limit t± «: t\\,f = 0; and (iii) the rung-fermion 
limit t x » fy, f = 0. In the paired limit, which we solved in 
Sec. IIVBI spinless fermions form correlated-hopping bound 
pairs, and so the ladder model can be mapped to the excluded 
chain of hardcore bosons. The ground state of this latter model 
was solved exactly in Sec.|ll]and its ground-state correlations 
calculated in Sec. [TTlJ By reinterpreting the excluded chain 
correlations in ladder terms, we realized that ladder SC cor- 
relations dominates at large distances over ladder CDW cor- 
relations, both of which decay as power laws with separation, 
with leading exponents 70 = \ and j8i = I + |(| - N2) re- 
spectively, N2 being the ladder density. We also showed, us- 
ing a restricted probabilities argument, that ladder FL corre- 
lations decay exponentially with separation, with a density- 
dependent correlation length. 

Next, in the two-leg limit, which we solved in Sec. IIVCI 
we argued based on a "particle-in-a-box" picture that succes- 
sive spinless fermions in the two-fold degenerate staggered 
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TABLE I: A summary of the leading correlation exponents and wave 
vectors of various correlation functions that decay as power laws in 
the (i) paired limit f » fy,/ ± ; (ii) two-leg limit t ± <<c /y, f = 0; 
and (iii) rung-fermion limit t ± » t\\, f = 0. The wave vector k of 
the leading terms in the correlation functions are reported in terms 
of k F = nN[ , where < N\ < \ is the excluded chain density. The 



suffixes n and 
model. 



indicate further symmetries possible in the ladder 
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FL 
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ground states occupy different legs of the ladder. We write 
these ground states exactly in terms of the one-dimensional 
Fermi sea in Sec. IIVC 11 before calculating correlations in 
Sec. IIVC 21 We found, using a different restricted probabili- 
ties argument, that FL correlations decay exponentially with 
separation, with a density-dependent correlation length. CDW 
and SC correlations symmetric (antisymmetric) with respect 
to a reflection about the ladder axis decay as power laws, with 
universal leading exponents B\ = 2(\) and j\ = 2(|) respec- 
tively. 

Finally, in the rung-fermion limit, we mapped the ladder 
model to an excluded chain of spinless fermions in Sec. lIVDl 
Since we have already solved this latter model in Sec. [II] and 
calculated its ground-state correlations in Sec. Hn]below half- 
filling (which corresponds to quarter-filling on the ladder), we 
discussed the phase separations that occurs on ladders with 
greater than quarter filling in Sec. II VP II Correlation expo- 
nents obtained for the three limiting cases of our ladder model 
Eq. dl.2| >, as well as those for the excluded chains of hardcore 
bosons and spinless fermions, are summarized in Table U 

In this study, we find the emergence of surprising univer- 



sal correlation exponents. In the Luttinger liquid paradigm, 
all correlation exponents can be written in terms of the 

. iQ in ii 

exponent o " " 



y p = \{K p + K p l -2), 
7cr = + K~ x - 2) 



(5.1) 



appearing in the quantum-mechanical propagator, also called 
the (equal-time) two-point function 



G(r) ~ Ai r a cosk F r = A\ r [1+2 ^' + ^ ,] cosk F r 
= A 1 r-H(W>+(W>] cos k F r. 



(5.2) 



The parameters K p and depend generically on the fill- 
ing fraction and the interaction strength, and thus all cor- 
relation exponents are non-universal. In particular, various 
theoretical approaches (see review by S61yom)i2 tell us that 
the charge density waves (CDW), spin density wave (SDW), 
singlet superconductivity (SSC) and triplet superconductivity 
(TSC) correlations decay as power laws 



<"(0)«(r)> 



+ B 2 r~ K "~ K ° cos 2k F r + 



B 4 r 4Kp cos4k F r, 



(o- x (0)o- x (r)) = (o- y (0)o- y (r)} 

5- + D 2 , xv r 

r z 

(cr z (0)o--(r)) ~^-+D 2 Z r- K "- K " cos 2k F r, (5.3c) 
r~ 



(5.3a) 

K "- K - cos2k F r, (5.3b) 



<Aj j0 (0)A<,,o(r)> = <A]" () (0)A 1 ,„(r)) ~ C r 

t 

i,±i 



(Aj ± ,(0)A li±1 (r)>~C r-^-^ 



(5.3d) 
(5.3e) 



in a Tomonaga-Luttinger liquid. 

When the chain of interacting spinfull fermions is spin- 
rotation invariant (for example, in the absence of an exter- 
nal magnetic field), the spin stiffness constant must be K a — 
1, and the ground-state properties become completely deter- 
mined by the single nontrivial Luttinger parameter K p . The 
spinfull power laws thus become 



G(r)~Air-K^ + *? i+ 2)cos*jpr, 

<n(0)n(r)> % + B-,r K ^ 1 cos2k F r + 

nr l 

B\r^ Kp cos4fc f r, 
1 



<o-(0) • air)) ~ —z + D 2 r 
nr z 



K »~ x cos 2k F r, 



<A>(0)A (r)> = <Al(0)A!(r)> ~ C r~^ 



(5.4a) 

(5.4b) 
(5.4c) 
(5.4d) 



For spinless fermions, there is only one independent stiffness 
constant K p - Kg- - so that the spinfull power laws 

which have proper spinless analogs are 



G(r) ~ Air - ^^^ 1 ) cos k F r, 

(n(0)n(r)) + B 2 r 2K cos 2k F r + 

nr z 



(5.5a) 
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B 4 r 4K cos 4k F r. (5.5b) 

In the Luttinger liquid paradigm, universal correlation ex- 
ponents only arise in the special case of the Fermi liquid, 
where we have K p = 1 . Consequently, the two-point function 
decays as G(r) ~ r _1 , while the CDW and SC correlations 
both decay as r~ 2 . However, the universal correlations expo- 
nents that we find in our exact solutions are different from 
these. Furthermore, the correlation exponents a, j3, and y of 
the FL, CDW, and SC correlations ought to obey definite re- 
lations in a Luttinger liquid, because they can all be written 
in terms of a single Luttinger parameter K. Again, the uni- 
versal and non-universal correlation exponents we find in our 
exact solutions do not obey these relations. These observa- 
tions bring us to the paper by Efetov and Larkin, who first cal- 
culated the universal FL correlation exponent for an ordinary 
chain of hardcore bosons to be a = \ M If we accept for the 
moment that the Luttinger paradigm is correct, and that uni- 
versal correlation exponents can only be found at the Fermi 
liquid fixed point, then we are led to the conclusion that a = | 
must be a correlation exponent of the Fermi liquid. Clearly, 
this exponent does not belong to the Fermi liquid FL corre- 
lation (which should be a = 1, so what correlation does it 
belong to? 

In the seminal paper by Jordan and Wigner, the or- 
dinary chain of hardcore bosons is mapped to the ordi- 
nary chain of spinless fermions using the Jordan-Wigner 
transformation (see Appendix [A}. In this transformation, 
the hardcore boson point operators b\ and bj +r are each 

mapped to spinless fermion string operators cj Y\i<j(~^) ni an d 

Y\i<j+ r (— l) n ' c 7+r respectively. The FL correlation (b^bj +r ) 
between two hardcore boson point operators thus become 
the expectation (cj tl'lj+i '( _1 )"'£/+/•) of the string operator 
ct n^j+i" (.~VF' c j+r- which Efetov and Larkin found to de- 
cay with separation r as r~'/ 2 . String correlations such as this 
have never been systematically studied. One reason for this 
lack of interest is that typical string correlations, which re- 
ceive contributions only from restricted classes of configura- 
tions, decay exponentially with r, as we have seen for the FL 
correlation in the paired limit (Sec. IIV B 2\ and the two-leg 
limit (Sec. lIVCZt . However, there appear to many string cor- 
relations that decay with separation r as power laws. These 
power law decays are associated with (quasi-)long-range or- 
der that we have not been creative enough to imagine. 

In Sec. IIV CT1 we found in the two-leg limit that the stag- 
gered ground state has long-range order, in that if we know 
the pth particle is on leg i = 1, then we know for certain that 
the (p + 2s)th particle is on leg i = 1, and the (p + 2s + l)th 
particle is on leg i = 2, even as s — > oo, and even though we 
have no idea where these particles are on the ladder. This 
long-range order is not the usual kind of long-range order, 
which can be written in terms of the correlation between local 
order parameters, but is a long-range string order. The map 
from the ordinary chain ground state to the staggered ladder 
ground state, which is the inverse of the one constructed in 
Sec. IIVC II implicitly involves string operators, in that if we 
take the pth particle in the ordinary ground state configura- 



tion, we will know whether to map it to a particle on leg i = 1 
or leg i = 2, after we know which legs the preceding particles 
are on. Also, while it is deceptively simple to describe what 
the string operator in this inverse map does, which is to project 
out any combination of more than or equal to two consecutive 
particles on the same leg of the ladder, we know of no com- 
pact way to write down the string operator, even in this simple 
limit, unlike for the case of the Jordan-Wigner string. 

What we do know, drawing parallels from the Jordan- 
Wigner map from hardcore bosons on ordinary chains to non- 
interacting spinless fermions, is that a string map from one 
model to another will map some products of local operators to 
string operators, for example, the hardcore boson b\bj +r to the 

spinless fermion ct n^lj+i 1 ( _ l)"' c /+/> an d other products of 
local operators to products of local operators, for example, the 
hardcore boson njrij+r to the spinless fermion « /«)+,-. Having 
understood this, we realized that the CDW+ and SC+ corre- 
lations in the staggered ground state get mapped to the the 
correlation of local operators, because the string operators in- 
volved in the map multiply and cancel each other. On the other 
hand, when we map the CDW- and SC- staggered ground- 
state correlations to correlations of a chain of noninteracting 
spinless fermions, the string operators involved in the map do 
not cancel each other, and thus the resulting ordinary chain 
spinless-fermion correlations are string correlations. We also 
realized that these string correlations are operationally defined 
by the intervening-particle expansions we used to compute 
them. 

Since all the exact solutions we have obtained in this pa- 
per can ultimately be mapped to the one-dimensional Fermi 
sea, we conjecture that all correlation exponents are univer- 
sal. We claim that: (i) all exponents that are explicitly uni- 
versal are simple rational polynomials of the single universal 
spinless Fermi liquid parameter K — 1 ; and (ii) non-universal 
exponents are the result of (under)fitting linear combinations 
of universal power laws to a single power law. For example, 
in the two-leg limit, the leading universal exponent y8i = | 
of the CDW- correlation in the staggered ground state can be 
shown using a bosonization calculation of the string correla- 
tion it is mapped to, to follow automatically from the universal 
Fermi liquid parameter K — 1 In this same limiting case, 
the leading universal correlation exponent y\ — % of the SC- 
correlation, which gets mapped to a significantly more com- 
plicated string correlation, can conceivably be written as the 
combination 



of the universal Fermi liquid parameter K = 1, even though 
the bosonized form of this string correlation is not known. For 
the excluded chain of hardcore bosons or spinless fermions, 
nonlinear curve fitting of the SC correlation to the sum of one 
leading power-law decay and one subleading power-law decay 
leads to weakly non-universal correlation exponents for both 
power laws, whereas a complicated sum of power-law decays, 
Eq. ( 13. 4K produces a better fit visually. We believe good fits 
can also be obtained, using similar complicated sum of power- 
law decays, for those numerical correlations which we found 
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to have strongly non-universal correlation exponents. 

Finally, we asked ourselves whether all these string corre- 
lations that we have predicted will decay with separation r 
slower than the two-point function (CjCj+ r ) ~ r _1 can be mea- 
sured in a chain of noninteracting spinless fermions. Since 
these string operators are nonlocal observables, they do not in 
general couple to local measurements, so direct experimental 
measurement would be challenging, if not downright impos- 
sible. However, we would like to suggest the following pos- 
sibility: for a given string correlation of the one-dimensional 
Fermi sea, cook up in the laboratory an experimental system 
in which the corresponding correlation is a point correlation. 
If the ground-state of the experimental system can be mapped 
to the one-dimensional Fermi sea, we expect a measurement 
of the point correlation exponent in the experimental system 
to be an indirect measurement of the string correlation expo- 
nent in the Fermi sea. 



APPENDIX A: JORDAN- WIGNER TRANSFORMATION 

On a one-dimensional chain, hardcore bosons cannot move 
past each other, as one boson must first hop on top of the other 
— a move explicitly forbidden by the hardcore condition — 
for this to happen. For a different reason (the Pauli Exclu- 
sion Principle), but to the same effect, noninteracting spinless 
fermions on a one-dimensional chain cannot exchange posi- 
tions. Therefore, in one dimension, the hardcore-boson and 
noninteracting-spinless-fermion Hamiltonians are also identi- 
cal in structure, and thus the ground state of a chain of hard- 
core bosons is related to the Fermi-sea ground state of a chain 
of noninteracting spinless fermions in a simple way. A trans- 
lation machinery exists to map back and forth between these 
two ground states. This is the Jordan-Wigner transformation 23 

bi = Y\{l-2c)cj)c h b\=c\\\{t-2c)c 3 ), (Al) 

j<i j<i 

which maps hardcore bosons to spinless fermions, where the 
product 

f]a - 2c]cj) = Y\(l - 2nj) = Y\(-l) nj (A2) 

j<i j<i ]<i 

is called the Jordan-Wigner string. 

In Section IIVBI we saw how pairs of spinless fermions 
bound by correlated hops in the limit t' » t±,t\\ can be 
mapped to hardcore bosons with infinite nearest-neighbor re- 
pulsion, and then to hardcore bosons using the right-exclusion 
map described in Section III Al and then finally to noninter- 
acting spinless fermions. In Section III Bl we saw how ex- 
cluded hardcore-boson expectations are related to appropri- 
ately chosen ordinary hardcore-boson expectations. This re- 
lation between excluded hardcore-boson expectations and or- 
dinary hardcore-boson expectations will typically involve the 
intervening-particle expansion Eq. ( 12.121 ). As such, we will 
encounter hardcore-boson expectations of the form 

{b\(t - ■ ■ ■ «,+/, ■ ■ ■ n i+ i p ••■(!- n i+ r>-\)bi +r >) , (A3) 



a lot, where there are p hardcore-boson occupation number 
operators «, + /, at sites i + I, and r' — p — 1 hardcore-boson 
operators (1 - «,+/'), at sites i + 1', between the hardcore boson 
operators b ' at site i and 

bi+r' at site i + r . 
To evaluate these expectations, we first invoke the Jordan- 
Wigner transformation (IAU to replace all the hardcore-boson 
occupation number operators tij = bjbj by spinless-fermion 
occupation number operators rij - c\cj in dA3l ). Then, to 
account for the two unpaired hardcore-boson operators at the 
ends of the hardcore-boson operator product, we write dA3l ) 
as the spinless-fermion expectation 

{c\ ]~ [(1 - In - «, + i) • • • n i+ i, ■ ■ ■ n i+lp •••(!- n,+,.<-i) x 
j<> 

f](l-2n / ) \\ (l-2nj)c i+ ,}. (A4) 

j<i i<j<i+r> 

Noting that all Jordan-Wigner string operators (1 - 2nj) com- 
mutes with riji and (1 - ny), for j < i and i < f < i + r', and 
that 

(1 - 2nj)(l - 2nj) = 1, (A5) 

we can bring the Jordan-Wigner string r[;</(l ~ 2n/) associ- 
ated with the annihilation operator c, +( - through the interven- 
ing spinless-fermion operators to obtain 

(cj(l-flj+i) ■ ■ • m+i x ■ ■ ■ hm p ■ ■ ■ (l-n !+r <_i) [ (t-2nj)Ci+r'). 

i< j<i+r' 

(A6) 

Then, using the fact that 

c + (l-2« ; ) = c + , nj(l-2nj) = -n h (t-nj)(t-2nj) = (!-«,), 

(A7) 

we can finally write the hardcore-boson expectation 

{b] Y](t- nj) Y] nj b l+ ,) = (-1)" {c] Y](t- nj) f] nj c i+r> ) 

empty filled empty filled 

(A8) 

as a spinless-fermion expectation, where p is the number of 
occupied sites between ; and i + r'. The suffixes 'empty' or 
'filled' in the products in (IA8t refer to the sites between i and 
i + r which are empty or filled respectively. 
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